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Abstract. In a cut-off Woods-Saxon (CWS) potential with realistic depth S-matrix poles being far from 
the imaginary wave number axis form a sequence where the distances of the consecutive resonances are 
inversely proportional with the cut-off radius value, which is an unphysical parameter. Other poles lying 
closer to the imaginary wave number axis might have trajectories with irregular shapes as the depth of the 
potential increases. Poles being close repel each other, and their repulsion is responsible for the changes of 
the directions of the corresponding trajectories. The repulsion might cause that certain resonances become 
antibound and later resonances again when they collide on the imaginary axis. The interaction is extremely 
sensitive to the cut-off radius value, which is an apparent handicap of the CWS potential. 

PACS. 21.10.Pc - 25.40.Dn 


1 Introduction 

Gamow shell model (GSM) [T] became a useful tool in 
analyzing drip line nuclei produced in laboratories with 
radioactive beam facilities. A most recent analysis of this 
type is that of the ^Be(p, 7 )®B and the ^Li(n, 7 )®Li re¬ 
actions in Ref. [5]. The key elements of GSM are the 
Berggren-ensembles of single particle states. The Berggren- 
ensemble might include resonant and sometime antibound 
states beside the bound states and the scattering states 
along a complex path. The shape of the path determines 
the set of the S'-matrix pole states to be included in the 
ensemble. In order to have smooth contribution from the 
scattering states the shape of the path should go reason¬ 
ably far from the poles. Therefore to know, where the poles 
are located has crucial importance. 

The most common phenomenological nuclear potential 
used is the cut-off Woods-Saxon (CWS) potential. Two 
handicaps of the CWS potential are that the positions of 
the broad resonances do depend on the cut-off radius and 
the pole trajectories also strongly depend on this radius, 
which is an unphysical parameter. Moreover some of the 
trajectories show strange circulating shapes [3]- Two of 
us suggested an alternative potential, the SV potential [3] 
which goes to zero smoothly at a finite distance. With 
SV potential the strange circulating shapes are missing, 
therefore it was conjectured [3] that the circulation is due 
to the jump of the CWS at the cut-off radius. In spite of 
of these handicaps the popularity of using CWS potential 
is still survives, therefore in this work we study further 


the pole distribution and the pole trajectories in CWS 
potential. 

2 Formalism 

Very little is known about the positions of the broad res¬ 
onances which correspond to poles lying farther from the 
real energy or wave-number axes. By varying the strength 
of the potential we can calculate pole trajectories on the 
complex wave number (fc) plane. The decaying resonances 
lie in the fourth quadrant of the k plane. The captur¬ 
ing resonances are mirror images of the decaying reso¬ 
nances in a real potential, therefore their trajectories are 
simply reflected to the imaginary fc-axis. The complex 
wave number of a Gamow resonance is k = + ik^ 

with k^ < 0, and k^ > 0 for a decaying resonance, and 
k^ < 0 for a capturing resonance. The energy is propor¬ 
tional to , therefore the unbound poles lie on the second 
energy-sheet. Resonances lying close to the imaginary k- 
axis, when —tt < arg[k) < —7r/2, i.e. |fc^| < |fc^|, are 
called as sub-threshold resonances [S]. We want to calcu¬ 
late poles of the S'-matrix in a given partial wave having 
orbital angular momentum /. The partial wave solution 
satisfies the radial Schroedinger equation 

u"{r, k) kf{r)u{r, k) = 0 , (1) 

where prime denotes the derivative with respect to (wrt) 
the radial distance r. The function 
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is the so called squared local wave number. 

In this work we consider a very simple case in which 
a spin-less neutral particle is scattered on a spherically 
symmetric target nucleus represented by a real CWS type 
nuclear potential v{r). The first boundary condition (BC) 
for the solution u{r,k) is its regularity at r = 0: 

m(0,/c)=0. (3) 

The other BC is specified at large distance in the asymp¬ 
totic region, at i?ass, i-e. beyond the cut-off radius Rmax, 
where 

v{Rass > Rmax) = 0 . (4) 

For a scattering state the asymptotic BC requires that the 
solution M(r, k) should be a combination of the incoming 
[kr) and outgoing H^{kr) free waves satisfying the 
Ricatti-Hankel differential equation: 

u(r, k) = A[H^ {kr) - S{k)H+ (fcr)] . (5) 

Here S{k) is the element of the scattering matrix in the 
Ath partial wave. 

For the CWS potential the radial equation in Eq. © 
can be solved only numerically. At r = Rass the numerical 
solution should match to that of the asymptotic equation 
in Eq. ®. We can calculate S{k) from the logarithmic 
derivative at Rass'- 


^i{Rass: k) 


^i{Rassi 
(.Rass ^ 


( 6 ) 


were Ui(r,k) is the internal solution, which is regular at 
r = 0. 

The value of the 5'-matrix at a given wave number k 
can be calculated as: 


(^^ass) ^i(R ass 1 k)H-(kR 

ass') 

kH+(kRass) - Z,(Rass,k)H+(kRass) ’ 
where bar denotes derivative wrt the argument kr. 


3 Resonances as poles of the 5'-matrix 


For a resonance the asymptotic BC is defined differently. 
It is required that the solution u{r, k) and its derivative 
u'(r, k) should be outgoing type and the logarithmic deriva¬ 
tive of the external solution is 


Ze(r, k) = k 


(kr) 

H+(kr) ■ 


( 8 ) 


Both BC (in Eqs. ([21) and (jS])) can be satisfied simultane¬ 
ously only at discrete complex k eigenvalues, at the poles 
of S(k). Here the complex k eigenvalue is fixed by the ze¬ 
ros of the difference of the logarithmic derivatives of the 
internal and the external solutions 


(9) 


The computer programs GAMOW [5], and ANTI [7] find 
the zeros of G(k,r), at certain Rm matching radius 0 < 
Rm < Rass- For a broad resonance the proper choice of 
this Rm is difficult. The zero is searched by Newton it¬ 
erations, and the iteration process often converges poorly 
or fails. Therefore we developed a new method in which 
we compare the logarithmic derivatives in a wide region 
in r. This method is built into the program JOZSCQ |S]. 
We calculate G(k, r) in Eq. ([2]) at equidistant mesh-points 
with mesh size h at rj = jh, j G [* 1 , 12 ]- The mesh points 
are taken from a region where the nuclear potential falls 
most rapidly. Then we search for the absolute minimum 
of the function of two real variables k^, and k^: 

f(fc«,fc') = logli;iG'(fc,ri)|| . (10) 

j=h 

Absolute minima of the function F(k^, k^) in Eq. (US 
should have a large negative value. The position of the 
absolute minimum is the pole position of S(k). The mini¬ 
mum of the function is found by using the Powell’s method 
in Ref. [2]. To find the minima of the function F(k^,k^) 
hrst we explore the landscape of the function F(k) in a 
complex k domain of our interest. We plot the function 
—F(k), which has peaks at the poles of S(k) on a grid 
(see e.g. Fig. |T|). 


4 Pole trajectories in CWS potential 

The nuclear potential in Eq. ([2]) is in fm“^ units, like the 
k'^ is. It is more usual however, to give it the same units in 
which the energy is given, i.e. in MeV. It is also usual to 
give the potential as a product of its strength Vb in MeV 
and its dimensionless radial shape: 

V^^^(r,R,a,Rma^) = -Ho/cws(r, R, a, R„,ax) , (11) 

where the later is 

/cWS (^; Ri Rmax) — ^(Rmax ^) r-R - (1^) 

1-be— 

Here 0(x) denotes the Heaviside step function, being zero 
for negative and unity for non-negative arguments. The 
parameters of the shape in Eq. ((T2]) are the radius R, the 
diffuseness a and the cut-off radius Rmax- 

Sometimes it is useful to draw the trajectories of the 
individual poles, i.e. to explore how the positions of the 
poles change when we keep the shape of the CWS poten¬ 
tial fixed but vary the Vq strength of the potential. By 
increasing Vb the pole trajectory moves to the upper half 
of the fc-plane, where it becomes a bound state with neg¬ 
ative real energy and with purely imaginary k. Here the 
solution becomes a real function with finite number of ze¬ 
ros. The nodes number n counts the zeros of the bound 
state wave function outside r = 0. 

^ The program name is chosen to honor the late Jozsef 
Zimanyi to whom one of the authors (T. Vertse) is grateful 
for starting his carrier. 


G(fc, r) = z^(r, k) - Ze(r, k) . 
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Fig. 1. Landscape of the —F{k) surface over the complex wave 
number domain € [—0.1, 2.5], and k^ £ [0.1, —0.7] for I = 2, 
and Umax — 15. 

Antibound solutions belong to poles of S{k) on the 
negative wing of the imaginary fc-axis. They diverge with¬ 
out oscillations in the asymptotic region as r —>■ oo. The 
poles of the S{k) lying off the imaginary fc-axis on the 
lower half of the fc-plane are resonances. For a resonance 
the wave function is complex and the tail of the wave 
function (in both the real and the imaginary parts) oscil¬ 
lates around r-axis with exponentially growing amplitude 
as r —>■ oo. Both parts have infinite number of zeros. 

It was shown earlier [3] that in the CWS potential the 
positions of broad resonances do depend on the value of 
I^max [HID], therefore, the cut-off radius is an important 
but unphysical parameter of the CWS form in Eq. 


5 Numerical example 

Here we study the positions of the poles for Z = 2, where 
we have a centrifugal barrier. Among the resonances we 
have poles with small imaginary parts. These narrow res¬ 
onances appear only if we have a barrier in the poten¬ 
tial. As example we calculate the neutron single-particle 
resonances in the ^°®Pb-|-n system in the CWS potential 
used in Ref. [12] but without the spin-orbit term. The 
values of the potential parameters are: Vb = 44.4 MeV, 
R = ro208i/3^ ro = 1.27, a = 0.7. Here we present only 
I = 2 results, which show a typical example for the I > 0 
cases. 

In Fig. [1] we present the —F{k) surface calculated for 
Rmax = 15 fm for CWS potential. The peaks of the surface 
show the approximate pole positions of the resonances. 
They give only a hint for the location of the poles, and 
help us to give reasonably good starting values for finding 
the poles of S{k) more accurately, when we calculate the 
minima of the quantity F{k) in Eq. (HOD for each pole. 
The positions of the three lowest lying poles in the domain 
0 < jfc^l < 0.8, and —0.7 < k^ < 0 are listed in Table [TJ 
Decaying resonances are denoted by dm with k^ > 0, and 


Table 1. Sequence numbers m, km values and node numbers 
n for the first three lowest lying poles for £ = 2, Rmax = 15 
fm, Vo = 44.4 MeV. 


m 

[fc^i 


n 

1 

0.076 

-0.147 

4 

2 

0.438 

-0.603 

3 

3 

0.752 

-0.449 

5 



Fig. 2. Trajectories of the first three poles as the depth of the 
CWS potential increases. Starting points of the trajectories at 
Vb = 44.4 MeV are denoted by circles. We use Rmax = 15 fm. 

capturing ones by Cm with k^ < 0- We assign a sequence 
number m to the poles as their jfc^j value increase. 

Note that the node number in the fourth column does 
not increase monotonously with m in the first column. 
Observe also that the poles with m = 1 and m = 2 are 
sub-threshold resonances. 

In the left corner of the landscape one can see a tiny 
pair of poles. They are the to = 1 poles in Table [T] (ci, di). 
On the right and middle part of the figure we can see a 
mountain with almost equidistant peaks. 

The ds resonance is a member of this peak-sequence 
lying closest to the imaginary axis. Below this d^ pole lies 
the d 2 sub-threshold resonance with considerably larger 
|fc^| value than the rest of the peaks forming the mountain. 

To understand the dynamics of the poles it is very use¬ 
ful to produce an animation (see Ref. |13j l in which the 
potential depth is changing by AVq in each temporal step 
and we can see how the landscape —F{k) changes in a 
given k domain as a function of Vq. In the animation we 
can conveniently follow the move of the poles as a func¬ 
tion of Vq. In Fig. [5] we show the trajectories of the poles 
with starting positions listed in Table [T] We plot their 
trajectories in a domain in the third and forth quadrants. 

By increasing Vb starting from Vb = 44.4 MeV the pair 
(ci,di) in Fig. [2] moves downward and the pair (c 2 ,d 2 ) 
upward. The third pair (csjda) moves toward the imagi¬ 
nary axis and interacts with the second pair ( 02 ,^ 2 ). The 
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Fig. 3. Repulsion of resonances d 2 and da as the depth of the 
CWS potential increases with steps AV = 0.1 MeV. Starting 
points of the trajectories at Vo ~ 44.4 MeV are denoted by 
circles. Their position is the closest at Vo = 51.5 MeV. 


repulsive interaction between them is shown in a magni¬ 
fied scale in Fig. [3] The ^2 and ds resonances repel each 
other when they come close at Vq = 51.5 MeV. The re¬ 
pulsion changes the direction of the ^2 trajectory slightly 
towards the imaginary axis, while ds turns down and move 
farther from the imaginary axis. The repulsion produces 
some sort of avoided crossing of the resonance levels. For 
narrow resonances similar phenomenon was observed in 
axially deformed Woods-Saxon potential in Ref. m when 
the deformation was changed. The repulsive interaction is 
extremely sensitive to the value of Rmax- A small change 
from Rmax = 15 fm to Rmax = 15.1 fm reduced the clos¬ 
est distance between the resonances to half of the distance 
shown in Fig. [S] 

Increasing the Vq value further the ( 02 ,^ 2 ) pair of poles 
interacts with the pair (ci,di) and pushes it to the imag¬ 
inary axis at about k = (0, —0.35). After the collision the 
pair {c 2 ,d 2 ) move upward and the C 2 pole collides in the 
origin with the ^2 pole. They form a double pole there. In¬ 
creasing the strength further the C 2 pole becomes to be a 
bound pole with node number n = 3. The other pole what 
we still call as d 2 moves down along the imaginary axis. 
(In identifying the separating poles from the double pole 
in the origin, we shifted them a bit from the origin by ap¬ 
plying a small imaginary potential strength.) At Vb = 66.4 
MeV d 2 antibound pole meets the upward moving Ci anti¬ 
bound pole at about k = (0, —0.1) on the imaginary axis. 
After the collision of the two antibound poles they move 
off the axis, ci turns left, c ?2 turns right and they become 
resonances again. The ci and ^2 resonances move along 
cucumber like paths symmetrically to the imaginary axis 
and at the ends of the paths they collide in the origin form¬ 
ing a double pole there. After that the ci pole becomes to 
be a bound state pole with node number n = 4, and ^2 
pole starts moving down again along the imaginary axis 
as an antibound pole. We shall see that the pole ^2 plays 


the role of a bouncer, since it creates bound state poles 
from the double-poles in the origin. 

An unexpected feature of the I = 2 case is that reso¬ 
nances meet with ^2 pole twice and they form double pole 
twice on the imaginary axis. First when the cucumber¬ 
like path starts and the second time when they meet in 
the origin. 

The same pattern is repeated again and again, i.e. the 
members of the symmetric pair of poles meet on the imag¬ 
inary axis well below the origin and the member coming 
from the forth quadrant turns down, the other coming 
from the third quadrant turns up and collides higher with 
the bouncer (^ 2 ) twice. 

6 Summary 

As a summary, we can conclude that for CWS potentials 
the I = 2 pole positions are sensitive to the value of the 
cut-off radius and their interaction is extremely sensitive 
to it. The distance of the poles in the mountain is regu¬ 
lated by the Rmax value, as it was observed in Ref. [3] for 
I = 5. The larger the Rmax value is, the closer the peaks 
are to each other. The close lying poles might interact 
with each other. When the resonances come close to the 
imaginary axis they might become antibound states and 
they might interact with other antibound poles moving 
along the imaginary axis. The interactions between poles 
distort the pole trajectories. Such distorted trajectories 
have been observed earlier in Ref. m where the start¬ 
ing values of the I = 0 trajectories were calculated for 
^®F. (Starting value of the trajectory is, when Vb is very 
small). In that paper it was observed that the positions of 
the starting values of the I — 0 resonances are governed 
by the Rmax value. In Fig. 8. of that reference the n = 1 
trajectory shows a loop. Strange circulating shapes were 
observed for I — 0 also in Fig. 1. of Ref. [3] for the n = 1 
pole trajectories of the CWS for different Rmax values. 

It is reasonable to assume that the repulsive interac¬ 
tion between poles is responsible for all these strange tra¬ 
jectory shapes. Large value of Rmax pushes the resonance 
pole too close to each other, where they interact. To re¬ 
duce the possibility of this interaction we suggest to use 
a relatively small value for the cut-off radius not much 
larger than the range of the best fit SV potential. 

Acknowledgement 

Authors are grateful to I. Hornyak for valuable discus¬ 
sions. This work was supported by the OTKA Grant No. 
K112962. 


References 

1 . N. Michel, W. Nazarewicz, M. Ploszajczak, and T. Vertse, 
J. Phys. G. 36, 013101 (2009). 

2. K. Fossez, N. Michel, M. Ploszajczak, Y. Jaganathen, and 
R. M. Id Betan, Phys. Rev. C 91, 034609 (2015). 





A. Baran et al.: Calculating broad neutron resonances in a cut-off Woods-Saxon potential 


5 


3. T. Vertse, R. G. Lovas, P. Salamon, and A. Racz, AIP 
Conf. Proc. 1491, 113 (2012). 

4. P. Salamon and T. Vertse, Phys. Rev. C 77 , 037302 (2008). 

5. A. M. Mukhamedzhanov, B. F. Irgaziev, V. Z. Goldberg, 
Yu. V. Orlov, and I. Qazi, Phys. Rev. C 81, 054314 (2010). 

6 . T. Vertse, K. F. Pal, Z. Balogh, Computer Physics Com¬ 
munications 27, 309 (1982). 

7. L.Gr. Ixaru, M. Rizea, T. Vertse, Computer Physics Com¬ 
munications 85, 217 (1995). 

8 . Cs. Noszaly, A. Baran, P. Salamon, and T. Vertse, JOZSO, 
a renewed version of the computer code GAMOW for cal¬ 
culating broad neutron resonances in phenomenological nu¬ 
clear potentials, to be submitted to Computer Physics 
Communications. 

9. W. H. Press, B. P. Flannery, S. A. Teukolsky, W. T. Vet- 
terling, Numerical Recipes, Cambridge University Press, 
Cambridge, 1987. 

10. J. Darai, A. Racz, P. Salamon, and R. G. Lovas, Phys. 
Rev. C 86, 014314 (2012). 

11. A. Racz, P. Salamon, and T. Vertse, Phys. Rev. C 84, 
037602 (2011). 

12. P. Curutchet, T. Vertse, and R. J. Liotta, Phys. Rev. C 
39, 1020 (1989). 

13. www.atomki.hu/atomki/TheorPhys/TrajAnimL2R15.avi 

14. L. S. Ferreira, E. Maglione, and R. J. Liotta, Phys. Rev. 
Lett. 78, 1640 (1997). 

15. P. Salamon, R. G. Lovas, R. M. Id Betan, T. Vertse, and 
L. Balkay, Phys. Rev. C 89, 054609 (2014). 


Published in Eur. Phys. J. A (2015) 51: 76. 
DOI 10.1140/epja/i2015-15076-l 




